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Discovery of VHE gamma-ray emission from the SNR G 15.4+0.1 with H.E.S.S. 
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Abstract: Supernova remnants (SNRs) have emerged as one of the largest source classes in very-high-energy (VHE; 
E>0.1TeV) astronomy. Many of the now known VHE 7-ray emitting SNRs have been discovered by the H.E.S.S. 
imaging Cherenkov telescope array, thanks to its unique access to the inner galaxy. Statistically-significant emission of 
VHE 7-rays has now been detected from the direction of the supernova remnant G 15.4+0.1. While the centroids of 
the H.E.S.S. source and the shell-type SNR are compatible, the VHE morphology suggests a center-dominated source at 
TeV energies, something which is at odds with the shell-like morphology observed at radio frequencies. This suggests 
that H.E.S.S. may be observing TeV emission from a previously unknown pulsar wind nebula (PWN) located within the 
boundaries of the radio shell. If this interpretation is correct, G 15.4+0.1 would in fact be a composite SNR, the first 
case in which an SNR is identified as a composite on the basis of VHE 7-ray observations. Archival data from MAGPIS 
gives exciting hints that there is radio emission from the central parts of the remnant, giving support to this hypothesis. 
Unfortunately, image artefacts from a nearby strong radio source produce considerable uncertainties in the radio analysis. 
Additional observations in both the radio and X-ray are needed to confirm the composite nature of G 15.4+0. 1 suggested 
by H.E.S.S. 
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1 Introduction 

Very-high-energy (VHE; E>0.1TeV) 7-ray emission has 
now been detected from more than 60 Galactic sources, 
and the vast majority of these are in some way connected 
to supernova remnants (SNRs) [1]. TeV 7-rays have been 
observed from SNR shells themselves and also from the 
pulsar wind nebulae (PWNe) at the center of many SNRs. 
While the former gives rise to VHE 7-ray emission tracing 
the shock-front of the SNR, the latter shows a centre-filled 
morphology, i.e. a central nebula, generally visible from 
low-energy radio wavelengths up to 7-ray photon energies 
of tens of TeV in the most extreme cases. Young PWNe 
typically have a rather symmetric morphology and are cen- 
trally located inside the host SNR. Such PWNe show a dis- 
tinct spectral fingerprint with a roughly flat radio spectral 
index and a much steeper spectrum at higher energies [2] . 
The pulsar powering the nebula is usually found close to its 
birthplace at the centre of the SNR in such systems. How- 
ever, only in a few cases has the pulsar actually been de- 
tected. One prominent example of a young PWN is found 
in the SNR G 0.9+0.1 which was first classified as a com- 
posite SNR on the basis of observations in radio [3] then 
subsequently detected in TeV 7-rays [4], before ultimately 



the energetic pulsar was found in a deep radio observation 
of the source [5], 

The SNR G 15.4+0.1 is a poorly studied object which 
was initially discovered in a 90-cm survey of the inner 
Galaxy conducted by the Very Large Array (VLA) [6]. 
The SNR was reported to have a shell-like morphology 
with a size about 14' x 15' and an average spectral in- 
dex a = —0.6 ± 0.2, indicating that the radio emission 
is dominated by non-thermal synchrotron emission from 
the shell. For the reasons outlined above, this source is 
a prime candidate for H.E.S.S. observations. Furthermore, 
since G 15.4+0.1 is larger than the H.E.S.S. point-spread- 
function (PSF; ~ 6'), morphological studies with H.E.S.S. 
have the capability to distinguish between a shell-type SNR 
or PWN-like origin of the VHE 7-ray emission, thus shed- 
ding more light on this little known SNR. 

2 H.E.S.S. observations and results 

The region around G 15.4+0. 1 was observed with H.E.S.S. 
between 2004 and 2010; the dataset consists primarily of 
observations from the H.E.S.S. Galactic Plane Survey [1] 
and offset observations of nearby sources. The dataset 
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in this region has a live-time of — 145h after standard 
H.E.S.S. quality selection [7], although the effective live- 
time at the target in question is significantly lower due to 
the large average offset of the pointings from the target. 
The dataset was analyzed using the standard H.E.S.S. anal- 
ysis [7] which employs the Hillas second moment method 
[8] for discriminating between 7-ray- and hadron-induced 
extensive air showers (EASs). Hard cuts, where at least 
200 photoelectrons are required in each recorded EAS im- 
age, were used instead of standard cuts (80 p.e.) in this 
analysis, in part because they give the improved angular 
resolution which is crucial for this analysis. To generate 
2D images, the ring background method [9] was used with 
an inner radius of 0.7°. All presented H.E.S.S. results have 
been cross-checked with an alternative analysis chain us- 
ing an independent calibration and 7-ray /hadron separation 
[10]. 

2.1 Detection and morphology 

Fig. 1 shows the 7-ray excess counts in the region around 
the SNR G 15.4+0.1. The image has been smoothed with 
a Gaussian kernel with standard deviation 0.06°. The sig- 
nificance of the source detection amounts to 7.0 u using 
an integration radius of 0.1°. The source centroid and ex- 
tension were determined by fitting a 2D Gaussian, con- 
volved with the H.E.S.S. PSF (0.076° for this analysis), 
to the unsmoothed 7-ray excess map. The best-fit cen- 
troid is at a J200 o = 18 h 18 m 3.4 s ± 4.6° tat and 5 J200 o = 
-15°27'54" ± 68^' tat (1 - 15.41°, b - 0.17°), and the 
source is thus given the identifier HESS J 1 8 1 8 — 154. In ad- 
dition to the statistics uncertainty in this fit, there is also a 
systematic uncertainty due to the pointing precision of the 
telescope array (20", [11]). Both the best-fit position and 
the best-fit extension (a = 0.071° ± 0.015°, 68% contain- 
ment radius), are indicated in Fig. 1. 

The results of this 7-ray morphological analysis are sum- 
marized in Table 1, together with updated values for 
the centroid position and average extension of the SNR 
G 15.4+0.1 in the radio. The latter have been estimated 
geometrically and the errors quoted represent the differ- 
ence between the results for the 20 cm image and the 
90 cm image. A comparison of these results shows that 
the centroids of G 15.4+0.1 and HESS J1818— 154 are 
compatible within l-a errors. Furthermore, the extension 
of HESS J1818 — 154 is significantly (at more than 2a) 
smaller than the average extension of G 15.4+0.1. 





HESSJ1818-154 


G 15.4+0.1 


aj20oo 

<5j2000 

Size 


18 h 18 m 3.4 s ± 4.6; tat 
-15°27'54"±68^' tat 
8.5' ± 1.8g tat 


18 h lg m 2 s ± 3s 

-15°28'32" ± 14" 
(14.4' x 9.6') ± 36" 



Table 1: Position and size of HESS J1818 — 154 and SNR 
G 15.4+0.1. 
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Figure 1: The VHE 7-ray excess (in units of 7-rays 
arcmin~ 2 ) around SNR G 15.4+0.1. The image has been 
corrected for the varying exposure across the field-of-view 
and smoothed with a 2D Gaussian with a width of 0.06°. 
The H.E.S.S. PSF for this analysis is shown in the bottom 
left inset. The color scale is chosen such that the blue-red 
transition occurs at roughly 4-cr significance. The best- 
fit centroid of the 7-ray excess is indicated by a cross, 
the size of which corresponds to the sum of both sta- 
tistical and systematic uncertainties. The best-fit exten- 
sion of HESS J 1 8 1 8 — 154 is represented as a circle. The 
green contours show the intensity (at 0.0175, 0.035 and 
0.0525 Jybeam -1 ) of the radio emission from the SNR 
shell, based on 90-cm VLA observations. 

2.2 Spectrum 

Spectral analysis was performed in a circular region cen- 
tered on the best-fit position of HES S J 1 8 1 8 - 1 54 and with 
a radius of 0.2°, enclosing —90% of the 7-ray emission. In 
this region, 116 excess 7-ray counts were found. The dif- 
ferential spectrum, shown in Fig. 2, is well fitted by a power 
law dN/dE = ^(E/TeV) -1 with a photon index T = 
—2.1 ± 0.3 s tat ± 0.2 syst and a flux normalisation at 1 TeV 
of 0o = (4.5±1.5 st at±0.9 syst )xl0- 13 cm- 2 s- 1 TeV- 1 . 
This corresponds to an integral flux F(> 1 TeV) = (4.0 ± 
1.6 8ta t ± O.Sgyst) x 10~ 13 cirr 2 s~\ equivalent to -1.8% 
of the flux of the Crab Nebula in the same energy range. 

3 G 15.4+0.1: A new composite SNR? 

Since the discovered 7-ray emission is found to be concen- 
trated in the centre of the SNR G 15.4+0.1 it is natural to 
consider HESS J1818 — 154 as a candidate TeV PWN asso- 
ciated with the SNR. In this scenario, the PWN may also 
be visible in the radio. Archival radio data from MAG- 
PIS [13] was used to investigate this hypothesis. Observa- 
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Figure 2: The top panel shows the differential energy spec- 
trum of HESS J1818 — 154. The error bars represent l-a 
statistical errors. The result of a power-law fit to the spec- 
tral data points is also shown. Events with energies be- 
tween 0.6 TeV and 13.5 TeV were used in the fit, and the 
bottom panel shows the residuals. 



tions of the region around G 15.4+0.1 are available at both 
20 cm and 90 cm wavelengths, where the PSF (here, half- 
power beamwidth) of the observations are 6.2" x 5.4" and 
24" x 18", respectively. The radio intensities observed in 
the field-of-view of G 15.4+0.1 at these two wavelengths 
are shown in Fig. 3. The SNR has a clear shell structure at 
both wavelengths, but notable differences in morphology, 
in particular in the central region. While the SNR centre 
appears largely void of radio emission at 90 cm, hints of a 
nebula with a diameter of ~5' are visible at 20 cm. This 
suggests there could be a varying spectral index across the 
SNR and putative PWN. A closer inspection of the 20-cm 
image reveals a strong source in the far north of the obser- 
vation that may have produced spurious emission and im- 
age artefacts. The radio analyses in this proceeding should 
therefore be considered preliminary. 

With these caveats noted, the integrated flux densities for 
both the central and shell components of SNR G 15.4+0.1 
were determined. The images were first convolved to a 
common beam size of 25" x 25" and point sources 1 , were 
removed by masking areas 1.5 times the beam size. The 
RMS noise was then derived from each (convolved) im- 
age and were found to be ergo = 20mJybeam~ 1 and 
(720 = 4mJybeam~ 1 , respectively. Integrating within the 
shell and central regions, denoted by the (white) ellipses 
and circle in Fig. 3, flux densities were measured and are 
given in Table 2. Uncertainties in the flux density include 
contributions from i) image noise, ii) calibration, and iii) 
wrong zero levels, as suggested in [14]. From the flux den- 
sities derived at 20- and 90-cm wavelengths, spectral indi- 



cies could be derived for the two regions; they are given in 
Table 2. The shell exhibits a steep spectrum (a ~ — 0.7), in 
agreement with [6] and indicative non-thermal synchrotron 
emission. However, the central region is found to exhibit 
a nearly flat spectrum (a ~ 0.3), supporting the central 
PWN hypothesis [2], and further adding to the case that 
G 15.4+0.1 may actually be a composite SNR. 



Region 


Flux (90 cm) 


Flux (20 cm) 


Index 


Centre 
Shell 


(0.8 + 0.2) Jy 
(14.0+ 1.0) Jy 


(1.2 + 0.2) Jy 
(4.6 + 0.5) Jy 


0.3 + 0.2 
-0.7 + 0.1 



Table 2: Derived flux densities and spectral indicies (S oc 
v a ) for the central and shell region of G 15.4+0.1. 



4 Discussion and conclusions 

Since a chance positional coincidence of the centroid of 
the VHE 7-ray emission with the center of the radio shell 
of G 15.4+0.1 seems highly unlikely, an association be- 
tween HESSJ1818-154 and the SNR is assumed here. 
The question is then whether the VHE emission originates 
from the shell of the SNR or from a putative PWN embed- 
ded inside the SNR. However, since the extension of the 
H.E.S.S. source is significantly smaller than the extension 
of the SNR as seen in the archival VLA data, it is clear 
that the bulk of the VHE emission originates from a yet 
unknown source located within the shell. The discovery 
of a radio or X-ray PWN counterpart to HESS J1818-154 
would confirm this scenario. High-quality X-ray observa- 
tions of the region are unfortunately lacking, but prelimi- 
nary analyses of 20-cm and 90-cm radio observations sug- 
gest that the SNR is composed of two distinct regions with 
different spectra. The central region is found to be charac- 
terised by a relatively flat radio spectrum, characteristic of 
the central PWN of a composite SNR [2], while the shell is 
characterised by a steeper spectrum, indicative of the non- 
thermal synchrotron emission expected from an SNR shell. 
G 15.4+0.1 thus appears to be an example of a composite 
SNR, similar to G 0.9+0.1, the first composite SNR to be 
detected at TeV energies. If the nature of G 15.4+0.1 can 
be confirmed, HES SJ1818 — 154 would be the first case in 
which a SNR is identified as a composite SNR on the basis 
of VHE 7-ray observations. 

Additional supporting evidence for a PWN origin of the 
TeV emission would be the detection of a pulsar within 
the central regions of the SNR. Unfortunately, no pulsar 
is currently known, in radio or X-rays, inside the shell 
of G 15.4+0.1, nor are there any known hard X-ray point 
sources. However, a typical young radio pulsar located at 
^9 kpc with an assumed luminosity L21 cm ~ 56 mJy kpc 2 
(i.e. similar to the median of high-luminosity, young, 
rotation-powered pulsars estimated by [12]) would have a 
flux density of S21 cm ~ 0.08 mJy, considerably lower than 



1. From SIMBAD, http://simbad.u-strasbg.fr/ 
simbad 
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Figure 3: Images of VLA 20-cm (left) and 90-cm (right) radio observations of G 15.4+0.1 in units of Jybeam -1 . The 
integration regions are indicated by white ellipses and a circle, respectively. The centroids of HESS J 1 8 1 8 — 154 and 
G 15.4+0.1 are indicated by a white and green cross. The extension of the cross corresponds to the difference in the 
estimated centroid positions for the 20-cm and 90-cm data. The images are shown in their original beam sizes, the size of 
which is indicated by the small green ellipses in the upper left corners. 



the sensitivity of the current VLA data. Beaming effects 
could also make the pulsar virtually undetectable in radio. 
If future observations of G 15.4+0.1 are able to detect the 
putative central pulsar associated with HESS J 1 8 1 8 — 154, 
it is expected to have a very high spin-down luminosity 
(E ~ 4 x 10 37 erg s, as seen in other TeV 7-ray emitting 
PWNe [5]. 

Little is known about the distance and age of G 15.4+0.1. 
Using the uncertain £ — D relationship [15], a distance of 
10 ± 3kpc is derived. Adopting lOkpc as the nominal dis- 
tance of the SNR would imply a Sedov age of ^9 kyr, as- 
suming an ambient ISM density of 1 cm" 3 . However, vary- 
ing the density by ±50% and taking into account uncertain- 
ties in the £ — D relationship, an age as low as 4 kyr can be 
derived. These age estimates would make SNR G 15.4+0. 1 
older than similar systems previously detected in VHE 7- 
rays. 
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